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Origin of Two Types of X-Ray Outbursts in Be/X-Ray Binaries. 

I. Accretion Scenarios 
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Abstract 

We propose the new scenario for X-ray outbursts in Bc/X-ray binaries that normal and giant outbursts 
are respectively caused by radiatively inefficient accretion flows (RIAFs) and Bondi-Hoyle-Lyttleton (BHL) 
accretion of the material transferred from the outermost part of a Be disk misaligned with the binary 
orbital plane. Based on simulated mass-transfer rates from misaligned Be disks, together with simplified 
accretion flow models, wc show that mass-accretion rates estimated from the luminosity of the normal 
X-ray outbursts are consistent with those obtained with advcction-dominatcd accretion flows, not with the 
standard, radiative-cooling dominated, accretion. Our RIAF scenario for normal X-ray outbursts resolves 
problems that have challenged the standard disk picture for these outbursts. When a misaligned Be disk 
crosses the orbit of the neutron star, e.g., by warping, the neutron star can capture a large amount of 
mass via BHL- type accretion during the disk transit event. We numerically show that such a process can 
reproduce the X-ray luminosity of giant X-ray outbursts. In the case of very high Be disk density, the 
accretion flow associated with the disk transit becomes supercritical, giving rise to the luminosity higher 
than the Eddington luminosity. 

Key words: stars: binaries, emission-line — Be, pulsars: individual (A 0535-1-262, 4U 0115+634), X- 
rays: binaries 



1. Introduction 

Bc/X-ray binaries are high- mass X-ray sources com- 
prised of a Be star and a neutron star in a wide and, in 
most cases, eccentric orbit. A Be star is a non-supergiant 
early-type star with a circumstellar disk formed by vis- 
cous diffusion of matter ejected from the star. This de- 
cretion disk is geometrically thin and Keplerian, as arc 
viscous accretion disks (Lcc et al. 1991), and is, in nor- 
mal condition, truncated by the tidal/resonant interaction 
with the neutron star (e.g., Okazaki & Negueruela 2001). 
Because of its eccentric orbit, the neutron star can cap- 
ture gas from the Be disk only for a short span of time 
during the close encounter with the Be disk. Such a brief 
interaction makes the Be/X-ray binaries transient X-ray 
sources. They show two types of X-ray outbursts: nor- 
mal (Type I) outbursts, of which the X-ray luminosity 
Lx ^ 10^^~^^erg s~^ and the interval is the orbital pe- 
riod, and giant (Type II) outbursts, which are significantly 
brighter (Lx > 10'^'^ erg s~^) and less frequent than normal 
outbursts (Stella et al. 1986; Negueruela et al. 1998). 

Thanks to long-term, multi-wavelength monitoring ob- 
servations of Bc/X-ray binaries, it is now widely accepted 
that normal outbursts arc triggered by the mass trans- 
fer from a tidally truncated Be disk at or near pcriastron 
passage and that giant outbursts are somehow associated 
with warping episodes of the Be disk. However, the details 



of the interaction between the Be disk and the neutron 
star and resultant accretion fiows remain elusive. 

In previous studies, the accretion disk which forms 
around the neutron star has been implicitly assumed to 
be modeled by the standard disk (Shakura & Sunyaev 
1973). It is geometrically thin and Keplerian, supported 
by rotation in the radial direction and by gas pressure in 
the vertical direction, in which the heating by viscosity 
is balanced with the cooling by radiation from the disk 
surface at each radius (e.g., Frank ct al. 2002; Kato ct al. 
2008). Several observed features, however, challenge the 
standard disk picture for Be/X-ray binaries. For instance, 
the spin- up episode of the X-ray pulsar in Be/X-ray bina- 
ries are, in general, seen only during X-ray outbursts (e.g., 
Bildstcn ct al. 1997 for A 05354-262 and GRO J1744-I-28; 
Parmar ct al. 1989; Wilson ct al. 2008 for EXO 2030+375), 
which often lasts for a much shorter time than the orbital 
period. This suggests that the accretion disk, if formed, 
is transient. Large and rapid flux variations in X-ray out- 
bursts also support the idea of transient accretion disks. 
In contrast, the standard disk theory suggests that the 
accretion disk is not transient but persistent, because the 
viscous- accretion timescale is generally much longer than 
the orbital period. 

It also remains unclear what mechanism triggers gi- 
ant X-ray outbursts. There is observational evidence for 
warping of the outer part of the Be disk before/during gi- 
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ant X-ray outbursts. This episode has most clearly been 
seen in 4U 0115-1-634 (Negueruela et al. 2001; Reig et al. 
2007), where the Ha emission line profile, which was usu- 
ally double-peaked, varied between a single peaked profile 
and a profile with absorption core, with the timescale of 
a year or so. Although single peaked profiles and those 
with absorption core are usually observed for Be stars seen 
from the polar and the equatorial directions, respectively, 
the spin axis of the Be star cannot vary with such a short 
timescale. Thus, this drastic line-profile variation indi- 
cates that the angle between the observer and the Be disk 
plane varied according to the precession of a warped Be 
disk (Negueruela et al. 2001; Reig et al. 2007). 

4U 0115-1-634 is not the only system that has shown 
the warping episode of the Be disk. Recent spectro- 
scopic monitoring observations of another Be/X-ray bi- 
nary A 0535-1-262 have detected very complicated changes 
in the Ha line profiles during and after the 2009 giant X- 
ray outburst (Moritani et al. 2011). Analyzing the vari- 
ability of several line profiles, Moritani et al. (2012) found 
that the observed features can be explained by a precess- 
ing, warped Be disk, which is misaligned with the binary 
orbital plane. 

We note that the association of the warping episode 
with giant X-ray outbursts favors systems with misaligned 
Be disks over those with Be disks coplanar with the or- 
bital plane. In misaligned systems, the warping can bend 
the disk outer part towards the orbital plane, enabling 
the neutron star to capture more material from the Be 
disk than it does without warping. On the other hand, in 
aligned systems where the Be disks are coplanar with the 
binary orbital plane, it is difficult to understand how the 
warping of initially coplanar disks can enhance the mass 
supply to the neutron star. Given that the eccentric orbit 
of Be/X-ray binaries is attributed to a supernova kick, the 
same mechanism is also likely to give rise to the misalign- 
ment between the Be star's spin axis and the orbital axis 
of the binary. Therefore, it is important to study dynam- 
ical interactions in misaligned systems, despite that there 
is little compelling evidence of the misalignment between 
the Be disk and the binary orbital plane in Be/X-ray bi- 
naries. 

In this paper, we examine accretion processes in Be/X- 
ray binaries, with particular attention to A 0535-1-262 and 
4U 01154-634, two best-studied systems with rather dif- 
ferent orbital parameters. Based on the study of these 
systems, we propose that normal and giant X-ray out- 
bursts in Bc/X-ray binaries are respectively caused by ra- 
diatively inefficient accretion flows (RIAFs) and Bondi- 
Hoyle-Lyttleton (BHL) accretion of the material trans- 
ferred from a Be disk misaligned with the binary orbital 
plane, and that giant outbursts occur when the tilted Be 
disk crosses the orbit of the neutron star in the direction 
of periastron. 

The paper is organized as follows: We first argue, in 
section 2, that the standard disk accretion picture does 
not work as a model for X-ray outbursts of Be/X-ray bi- 
naries. Then, in section 3, we show that RIAFs naturally 
explain characteristic features of normal X-ray outbursts. 



In section 4, we point out that the accretion rate and 
timescale evaluated on the basis of our scenario for giant 
X-ray outbursts arc consistent with those observed. We 
also discuss in this section the possibility of supercritical 
accretion flows onto the neutron star in Be/X-ray binaries. 
Section 5 is devoted for discussion about observational im- 
plications of our scenarios, and section 6 concludes the 
paper. 

2. Classical Picture for X-ray Outbursts 

As summarized in the next section, the normal X-ray 
outbursts last only for a small fraction of the orbital pe- 
riod. This short timescale challenges the classical, stan- 
dard disk picture for Be/X-ray binaries. In this section, 
we first show that the standard disk model has too long 
accretion timescale to explain the duration time of nor- 
mal X-ray outbursts. Afterward, we discuss a couple of 
issues, where the standard disk picture is inconsistent with 
observations. We present alternative scenarios for X-ray 
outbursts in later sections. 

2.1. Viscous Timescale 

In the standard disk, accretion occurs in the viscous 
tinicscale t^is, which is given by 

(e.g., Frank et al. 2002; Kato et al. 2008), where a is 
the Shakura-Sunyaev's viscosity parameter, r is the radius 
measured from the neutron star, H is the vertical scale- 
height of the accretion disk at r, and ilx = {GMx /r^Y^'^ 
with AIx being the mass of the neutron star. A steady 
a-disk solution has the scale-height given by 

= 5.6 X lQ-^a,\"'m%^'m~'J'f^'\ (2) 

where ao.i = a/0.1, rhie — M/10^^gs~^, mi. 4 ^ 
Mx/1.4M0, and f is the radius normalized by the 
Schwarzschild radius rg = 2GMx/c? « 4.1 x lO^mi.4 cm. 
Substituting equation (2) into equation (1), we have 

ivis = &.2aQ\^^m^Q^^'^m\{lf^/'^ s. (3) 

In order to evaluate equation (3), we need to estimate 
the outer radius of the accretion disk and the mass- 
accretion rate there. In the next subsection, we esti- 
mate tvis by utilizing results from numerical simulations 
of tidal interaction in the Be/X-ray binaries A 0535-1-262 
and 4U 0115-1-634. 

2.2. Accretion Timescale Estimated from Numerical 
Simulations 

2.2.1. Numerical Method 

The simulations presented below were performed with a 
3D SPH code, which is a particle method that divides the 
fluid into a set of discrete "fluid elements" (=particles), 
and is flexible in setting various initial configurations. The 
code is the same as that used by Okazaki et al. (2002) (see 
also Bate et al. 1995). Using a variable smoothing length. 
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Orbital period Porb (d) 


110.2 « 


24.3 5 


Orbital eccentricity e 


0.47 9 


0.34 5 


Mass ratio q 


0.056 


0.074 


Semi-major axis a (cm) 


2.0 X 10^3 


6.7 X 10^2 



1: Giangrande et al. (1980), 2: Negueruela & Okazaki (2001), 3: Vacca et al. (1996), 4: Coe et al. (1975), 5: 
Rappaport et al. (1978), 6: Makishima et al. (1999), 7: Ferrigno et al. (2009), 8: Moritani et al. (2010), 9: 
Finger et al. (1994) 



the SPH equations with a standard cubic-spline kernel are 
integrated with an individual time step for each particle. 

In our code, the Be disk is modeled by an ensemble of 
gas particles with negligible self-gravity. For simplicity, 
the gas particles are assumed to be isothermal at 0.6 Toff 
with Teff being the effective temperature of the Be star 
(Carciofi & Bjorkman 2006). We model the mass ejec- 
tion process from the star by injecting gas particles at 
a constant rate just outside the equatorial surface of the 
star. Once injected, gas particles interact with each other. 
As a result, most of the injected particles fall back on to 
the Be star by losing the angular momentum and a small 
fraction of particles drift outwards, obtaining the angular 
momentum from the other particles. 

On the other hand, the Be star and the neutron star are 
represented by sink particles with the appropriate grav- 
itational mass. Gas particles that fall within a specified 
accretion radius are accreted by the sink particle. We 
assume that the accretion radius of the Be star is equal 
to the stellar radius while for the neutron star, we 
adopt a variable accretion radius of tl, where tl is the 
Roche-lobe radius given approximately by 

0.49g2/3 

n = D (4) 

^ 0.69g2/3 + ln(l + 5i/3) ^ 

(Eggleton 1983) with the mass ratio q = Mx/M^, where 
is the mass of the Be star, and the instantaneous dis- 
tance between the stars, D. 

In simulations shown in this paper, we adopt artificial 
viscosity parameters agpH = 1 or /3sph = Oj which roughly 
corresponds to the Shakura-Sunyaev viscosity parameter 
a = 0.1 except in very-low density regions (e.g., Okazaki 
et al. 2002). We set the binary orbit on the x-y plane with 



the major axis along the x-axis (the apastron is in the +x- 
direction). At t = 0, the neutron star is at apastron. It 
orbits about the Be star primary with the orbital period 
Porb and the orbital eccentricity e. 
2.2.2. A 0535+262 

A 0535-1-262 consists of an 09.7IIIe star and an X-ray 
pulsar with a relatively wide [Porb = 110. 2d (Moritani et 
al. 2010)] and eccentric [e = 0.47 (Finger ct al. 1994)] 
orbit. As the mass and radius of the Be star, we take 
25 Mq and 15 i?© from Vacca et al. (1996). With these 
parameters, the semi-major axis is a = 2.0 x 10^^ cm. The 
magnetic field strength of the neutron star is 4.3 x 10^^ G 
from the observed cyclotron resonance scattering features 
(Makishima et al. 1999). Table 1 summarizes system pa- 
rameters used in this paper. 

We have performed two simulations for A 0535+262. In 
one simulation, the Be disk is coplanar with the binary or- 
bital plane, while in the other simulation, it is misaligned 
with the orbital plane by 45° about the semi- minor axis. 
Each simulation has run until the Be disk is fully devel- 
oped and start showing a regular orbital modulation in 
the mass capture rate by the neutron star. The num- 
ber of gas particles at the end of coplanar and misaligned 
disk simulations are - 170,000 (t^SOPorh) and ~210,000 
{t = 50 Porb), respectively. 

Figure 1 shows the result from (a) the coplanar disk 
simulation and (b) the misaligned disk simulation. In the 
upper panels, the thick (blue) lines denote the mass cap- 
ture rate by the neutron star as a function of orbital phase. 
We note that in the coplanar case, mass capture rate has 
a single peak at the periastron, while it is multi-peaked 
when the Be disk is highly inclined from the binary orbital 
plane. This is because in misaligned disks, the neutron 
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star makes close disk encounter twice an orbit at orbital 
phases different from the periastron. We also note that 
the neutron star captures more mass from coplanar Be 
disks than from misaligned ones. This is because the dis- 
tance between the neutron star and the Be disk in the 
coplanar system is shorter than that for the misaligned 
Be disk rotated around the binary semi-minor axis. 

In the coplanar disk simulation, the peak mass-capture 
rate is M « 10^^p_ii g s~^, while in the misaligned simu- 
lation it is M PS 2 — 3 X 10^^/9_iig s~-^, where is the 
base density of the Be disk normalized by 10~^^gcm~^, 
a typical value for Be stars. As the right axis shows, 
these rates would produce the peak X-ray luminosity 
of ~ 2 X 10^^p_ii, erg s"-'^ in the coplanar case, and 4 x 
lO'^^yO-ii erg s~^ in the misaligned case, if the accretion 
time of the captured material on to the neutron star were 
negligible. Here, we have calculated the X-ray luminos- 
ity by Lx = GMxM/Rx « 1-0 x lO-^migiedd with Ledd 
being the Eddington luminosity given by 

iiTcGMx 



-^Edd — 



1.8 X ergs" 



(5) 



is the opacity of the 
= 10^ cm as the radius of 



where c is the speed of light and 
electron scattering. Wc take Rx 
the neutron star. 

Since the Roche-lobe radius tl is larger than the neu- 
tron star radius by many orders, the accretion time from 
this radius to the neutron star surface could significantly 
reduce the peak amplitude and broaden the shape of the 
accretion rate. It is thus interesting to estimate the accre- 
tion rate at the neutron star radius. For this purpose, 
wc first calculated the circularization radius Tcirc from 
the specific angular momentum of captured material, j, 
by assuming that rcirc is approximately given by rcirc = 
j^/(GMx). We then evaluated the viscous timescale tvis 
at r = Tcirc, using equation (3). In our simulations, the 
mass- weighted average of rcirc is 
(coplanar case) and rdrc 



3xl0"cm~0.02a 



7 X 10^° cm - 0.004a (mis- 
aligned case). The lower panels of figure 1 shows the 
resultant tvis as a function of orbital phase. In both simu- 
lations, the viscous (i.e., accretion) timescale turns out to 
be longer than the orbital period at most phases, in par- 
ticular, <vis 5 — 6Porb at the peak mass-transfer phase. 

Using the phase dependence of mass capture rate and 
viscous timescale, we estimated the accretion rate on to 
the neutron star as follows. With the mass capture rate 
Mciip{t) at t, we approximate the amount of mass captured 
for a short span of time between t = to and t = to + At by 
Mciip{to)At. As a very rough approximation of the viscous 
accretion process, we assume that this amount of mass is 
accreted by the neutron star over ivis between to + ivis / 2 
and to + 3tvis/2 at a constant rate of Mcap(to)^V^vis- 
Note that in this treatment the accretion timescale is 
taken as being independent of the mass capture/accretion 
history. The method is thus expected to provide a rea- 
sonable estimate of the accretion rate, when the accretion 
timescale increases with orbital phase or slowly decreases 
with phase with the changing timescale longer than the 
accretion timescale itself. However, it might create an 
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Fig. 1. Simulation result for A 0535-1-262 with (a) a Be 
disk coplanar with the binary orbital plane and (b) a Be 
disk misaligned with the orbital plane by 45° about the 
semi-minor axis. Upper panels: the orbital modulation in 
the mass capture rate by the neutron star [thick (blue) lines] 
and estimated mass-accretion rate at the neutron star ra- 
dius [thin (red) lines]. The right axis shows the X-ray lu- 
minosity corresponding to the mass-accretion rate. Lower- 
panels: the viscous timescale at the circularization radius 
of captured material. The data is folded on the orbital 
period over 25Porb £i * 5; SOPorb in the coplanar case and 
45Porb ^ i < SOPorb in the misaligned case. The mass-cap- 
ture/accretion rate and the X-ray luminosity are measured 
in units of 10^®p_iigs^^ and 10^^p_ii erg s~^ , respectively, 



-11 is the base density of the Be disk normalized by 

-3 



where p_ 

10~ii gem""*, a typical value for Be stars. The dashed line 
in each lower panel denotes the mass- weighted average of the 
viscous timescale. 



artificial, spiky feature in the resulting accretion rate in 
case the accretion timescale rapidly decreases with or- 
bital phase, where in real systems the rapid accretion of 
newly captured material will certainly be prevented by 
an already-existing accretion disk with a longer accretion 
timescale. 

Applying the above procedure to the phase-dependent 
mass capture rate and viscous timescale, we obtained the 
orbital modulation of accretion rate estimated at the neu- 
tron star radius. The result is shown by the thin (red) lines 
in the upper panels of figure 1. As expected from the fact 
that the viscous timescale is much longer than the orbital 
period, the obtained mass accretion rate is mostly con- 
stant and exhibits only short and small enhancements at 
phases with short accretion timescales. These simulated 
features contradict with the observational characteristics 
of normal X-ray outbursts, which exhibit a large and rapid 
change in the X-ray flux. 
2.2.3. 4U 0115+634 

4U 0115-1-634 consists of a B0.2Ve star and an X-ray 
pulsar with a relatively narrow (Porb = 24.3 d) and eccen- 
tric (e = 0.34) orbit (Rappaport et al. 1978). As the mass 
and radius of the Be star, wc take 19 and 8Rq from 
Vacca et al. (1996). With these parameters, the semi- 
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(a) 4U 0115+634 with a (b) 4U 0115+634 with a 
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Fig. 2. Simulation result for 4U 0115+634 with (a) a copla- 
nar Be disk and (b) a Be disk misaligned with the orbital plane 
by 45° about the semi-minor axis. The data is folded on the 
orbital period over 45Porb ^ SO^orb fof both simulations. 
The format of the figure is the same as that of figure 1. 

major axis a is a = 6.7 x 10"'^^ cm, which is about three 
times smaller than that of A 0535+262. These system 
parameters are summarized in Table 1. 

As for A 0535+262, we have performed two simulations 
for 4U 0115+634, i.e., a coplanar disk simulation and a 
simulation with a Be disk misaligned by 45° about the 
semi-minor axis. Both simulations have run over 50Porb 
until the Be disk is fully developed and start showing a 
regular orbital modulation in the mass capture rate by 
the neutron star. The number of gas particles at the end 
of coplanar and misaligned simulations are ~ 140,000 and 
~ 110,000, respectively. 

Figure 2 shows the result from (a) the coplanar disk 
simulation and (b) the misaligned disk simulation. The 
format of the figure is the same as that of figure 1. From 
the upper panels, we note that the peak mass-capture rate 
in this system is several times lower than in A 0535+262. 
This is because the tidal truncation of the Be disk is 
more efficient in shorter orbital period and smaller or- 
bital eccentricity (Okazaki & Negueruela 2001). Owing 
to a larger gap between the neutron star orbit and the 
Be disk, the neutron star of 4U 0115+634 has more dif- 
ficulty in capturing gas from the Be disk than the neu- 
tron star of A 0535+262 does. By the same reason, the 
pre-periastron peak, which is caused by the pre-periastron 
close encounter of the neutron star to the Be disk, disap- 
pears in the misaligned simulation of 4U 0115+634. The 
stronger tidal interaction in 4U 0115+634 also gives rise 
to longer viscous timescale (lower panels), because the 
captured material has larger specific angular momentum. 
As a result, in 4U 0115+634 the mass-accretion rate esti- 
mated at the neutron star surface (red lines) are much 
lower and less variable than in A 0535+262 (the prc- 
periastron spike in Fig. 2a is an artificial feature arising 



from our simplified treatment). Apparently, these simu- 
lated features are inconsistent with the observational char- 
acteristics of normal X-ray outbursts. 



2.3. Direct Accretion vs. Propeller Regimes 

Although we have shown in the above sections that the 
standard disk scenario does not agree with the observed 
short duration and large and rapid X-ray fiux variations 
of normal X-ray outbursts, one may argue that the onset 
of propeller regime can truncate the accretion regime and 
hence make the standard disk model still viable. Given 
a strong positive correlation between the spin period of 
the pulsar and the orbital period of the binary seen in 
Be/X-ray binaries (e.g., Corbet 1986), it is likely that the 
propeller effect of the rotating magnetosphere of the pul- 
sar plays an important role in systems with short orbital 
periods. Indeed, in 4U 0115+634, Campana et al. (2001) 
witnessed a very rapid increase of the X-ray flux to be at- 
tributed to the transition from the propeller regime to the 
direct accretion regime. It is not clear, however, whether 
the same mechanism plays an important role in the X- 
ray behavior of systems with much longer orbital periods, 
such as A 0535+262. 

Below we examine whether the propeller effect also 
works for A 0535+262, truncating the accretion regime 
short and making a large difference between the quiescent 
and outburst X-ray fluxes. For this purpose, we compare 
the corotation radius, re, where the local Keplerian fre- 
quency is equal to the rotation frequency of the magneto- 
sphere, with the magnetospheric radius, tm, obtained by 
equating the magnetic field pressure of the neutron star 
to the gas pressure of the accretion disk. For tm < fc the 
velocity of rigidly rotating magnetic field is sub-Keplcrian 
at the magnetospheric radius, which permits the mate- 
rial in the accretion disk to fall on to the neutron star 
along magnetic field lines. The system is then in the di- 
rect accretion regime. On the other hand, for tm > rc, 
the magnetosphere rotating at a super-Keplerian speed 
expels disk material, prohibiting the accretion to occur. 
This regime is called the propeller regime. 

The corotation radius rc is given by 



47r2 



1/3 



1.4xlO«Fi;/>;/f cm 



(6) 



here Pspin is the spin period of the neutron star. 



Assuming the magnetic field of the neutron star to be a 
dipole-like field, the magnetospheric radius tm is written 

as 

tm ~ 4.9 X IQ^km^^^'' m^]/'^ cm (7) 

(e.g., Frank et al. 2002). Here, /130 = /x/lO^o [Gcm^] is 
the normalized magnetic moment and is a constant of 
order unity which depends on the physics and geometry 
of accretion. For spherical accretion, k ^ 1. For disk ac- 
cretion, however, the previous results range from k ^ 0.5 
to A: ~ 1. In this paper, we adopt fc = 1 as a conservative 
value for disk accretion, which gives rise to an upper limit 
of magnetospheric radius for given m. The system is in 
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Fig. 3. The long-term variation of the ratio of the mag- 
netospheric radius tm to the corotation radius tq for 
A 0535+262. The magnetospheric radius is estimated from 
the observed fiux with Swift with a conservative value of fc = 1 . 

the accretion regime when the normahzed accretion rate 
satisfies the fohowing condition: 



mi6 > 1.4 X lO^fc'/'n^ifmii/ Vic- 



es) 

For disk accretion in A 0535+262, where Pspin = 103s and 
M30 = 4.3, we have rhig ^ 0.5 for direct accretion regime. 

In figure 3, we compare tm with using the mass 
accretion rate estimated from the observed X-ray flux 
of A 0535+262. We utilized the X-ray observations by 
Swift/BAT (15 - 50 keV) ^, converting its count rates to 
the energy fluxes. We assumed that 1 Crab = 1.386 x 10~* 
ergcm~^s~^. modeling the Crab spectra with nn — 
0.26 X 10^^ cm~^, r = 2.1, and power-law normalization 
being 10.0. Figure 3 shows the ratio of the magnetospheric 
radius tm to the corotation radius rc- The figure clearly 
shows that the system is almost always in the direct accre- 
tion regime, except for periods of very low X-ray fluxes. 
Therefore, we conclude that in A 0535+262, the large and 
rapid flux changes of normal X-ray outbursts with short 
duration are due to an intrinsic nature of the accretion 
flow, not controlled by the propeller effect. This safely 
rules out the standard disk scenario for normal X-ray out- 
bursts. 

3. Origin of Normal X-ray Outbursts 

The observed features of normal X-ray outbursts are 
summarized as follows: They occur at or near the peri- 
astron passage. The X-ray luminosity increases by about 
one order of magnitude with respect to the pre-outburst 
state {Lx ^ 10'^^"'^^ erg s~^). The duration time is a small 
fraction of the orbital period, typically 0.2 - 0.3 Porh- 
Sometimes the precursors and flares occur prior to the 
main outbursts (Klochkov et al. 2010a and Klochkov et 
al. 2010b for EXO 2030+375; Caballero et al. 2008 for 
A 0535+262; Campana et al. 2001 for 4U 0115+634). 
Among those, the flares are interpreted to be caused 



by the instability of magnetosphere (e.g., Postnov et al. 
2008). 

Although it is widely accepted that normal X-ray out- 
bursts are triggered by the mass transfer from the Be disk 
at or near periastron, the accretion mechanism still re- 
mains an open question. As shown in the previous sec- 
tion, the standard disk scenario fails to explain the obser- 
vational features noted above. Below we propose that the 
normal X-ray outbursts are caused by RIAFs on to the 
neutron star from a tidally truncated Be disk. 

3.1. Radiatively Inefficient Accretion Flows 



In the context of the standard disk model, the viscous 
heating rate per unit area, Q^, is assumed to balance with 
the radiative cooling rate per unit area, Qrad, in the energy 
equation. This approach is justified when the accretion 
disk is geometrically thin so that Qadv ^ {H/r^Qj^ can 
be neglected, where Qadv represents the advective cooling 
rate per unit area. However, in hot accretion flows with 
H/rKi Cs/vK ~ 1, where Cs and vk are the sound speed and 
Keplerian velocity of the disk, respectively, the advection 
cooling cannot be neglected in the energy equation. 

Abramowicz et al. (1995) found that there are two fam- 
ilies of advection dominated accretion flow solutions. One 
family consists of optically thin, advection-dominated ac- 
cretion flows (ADAFs), while the other consists of opti- 
cally thick, advection-dominated accretion flows, so-called 
slim disks. Both families are thermally and viscously sta- 
ble. From figure 3 of Abramowicz et al. (1995), ADAF 
solutions exist if to = A//(LEdd/c^) ^ 1 (or equivalently 
^ 2OTO1.4) whereas there are slim disk solutions if 
m ^ 10 (or rhie ^ 2OOTO1.4). In the following, we con- 
sider ADAF solutions as the RIAF model for normal X- 
ray outbursts. We will also discuss slim disk solutions in 
Section 4.4. 

3.1.1. Self-similar ADAFs Solutions 

Narayan & Yi (1994) derived a set of solutions of self- 
similar, optically thin ADAFs with Qadv = fQ+, where / 
exhibits the extent to which the flow is advection dom- 
inated. These flows give us most of important proper- 
ties of the general solutions, although Qiad is not ex- 
plicitly shown in their solutions. Therefore, we model 
accretion flows onto the neutron star by the self-similar 
ADAFs when they have low accretion rates to ^ 1 (or 
rhiQ ^ 2OTO1.4). By adopting the self-similar scaling law 
proposed by Narayan & Yi (1994) and Narayan & Yi 
(1995), we obtain a set of solutions for a ^ 1 as 

f ' cm s , 



« -2.1 X 10 ciao.i^ 
if w 4.2 X 10^C3^^mi.4r cm. 



2.1 X 10" 



-1/2^ 



rw2.7x lO^^cgf^^K, 



(9) 
(10) 

ll6aol™^4^~^^^gcm"^ (11) 

(12) 



http:/ /swift. gsfc.nasa.gov/docs/swift/rcsults/transicnts/ 
lA0535p262/ 



where Wr is the radial velocity, p is the density, T is the 
disk temperature, and ci « 0.53, C2 ~ 0.34, and C3 w 0.35 
are numerical constants. 
3.1.2. Observational Implications 

The infall timescale of the ADAF is evaluated by 
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coplanar Be disk 



(b) A 0535+262 with a mis- 
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Fig. 4. ADAF model for A 0535+262 with (a) a coplanar Be 
disk and (b) a Be disk misaligned with the orbital plane by 
45° about the semi-minor axis. The format of the figure is 
the same as that of figure 1, except that the thin (red) lines 
in the upper panels denote the estimated mass-accretion rate 
at the neutron star radius using the ADAF solution and that 
the lower panels show the infall timescale of the ADAF at the 
circularization radius of captured material (solid lines) and 
the mass-weighted average of the ADAF timescale (dashed 
lines). 
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Fig. 5. ADAF model for 4U 0115+634 with (a) a coplanar 
Be disk and (b) a Be disk misaligned with the orbital plane 
by 45° about the semi-minor axis. The format of the figure 
is the same as that of figure 4. 



^ADAF — 1 

\Vr\ 



' 2.0 X IQ-^miAC^^aQ lf^/^ s. 



(13) 



Applying equation (13) to the same simulation data as in 
figures 1 and 2 for A 0535+262 and 4U 0115+634, respec- 
tively, we have calculated the mass-accretion rate profiles 
in the case of ADAFs. The results arc shown in figure 4 
for A 0535+262 and figure 5 for 4U 0115+634. From the 



results for A 0535+262, we note that the ADAF solutions 
provide X-ray profiles with similar duration and amplitude 
to those of observed normal outbursts. The pre-periastron 
spike of the mass-accretion rate in the coplanar simulation 
is also reminiscent of an observed short spike prior to peri- 
astron (e.g., Camero-Arranz et al. 2012). Comparing the 
mass-accretion rate in the upper panel with the iADAF 
plot in the lower panel, we note that the spiky feature 
is caused by the material with very low specific angular 
momentum. 

The simulated accretion rate for 4U 0115+634 also has 
the accretion timescale in agreement with that of the ob- 
served normal outbursts. The X-ray luminosity obtained 
for the typical Be disk density (10~^^gcm~'^ at the base 
of the disk), however, seems too low to be consistent with 
the observations. This might explain why this system has 
shown normal X-ray outbursts only when the Be disk was 
fully developed. 

Next, we evaluate the Bolometric X-ray luminosity of 
ADAFs. For simplicity, we assume that the X-ray ra- 
diation emitted from the ADAF is produced only by 
the Bremsstrahlung cooling. The cooling due to the 
other processes such as the Synchrotron radiation and 
Inverse Compton process will be discussed in a forthcom- 
ing paper in the context of high energy emissions. The 
Bremsstrahlung cooling rate per unit surface area is given 
by 

Qbr= ebrp'T^/^iJ 

« 7.7 X lO^^q^ag^TOj^jmier^^/^ergs^^cm"^, (14) 

where ebr ~ 1-2 x 10^^ergs~^cm~^ is the emissivity of the 
Bremsstrahlung radiation. Then, the Bolometric X-ray 
luminosity is given by 



^ADAF = 27rrs" 



1.7 X lO^^ag ifc"i/2mi;4 fi^^ ' m 



fQhi-df 

4/7 -2/7 ■ 8/7 



ergs 



(15) 



for tm <. ' circ ■ 

On the other hand, the X-ray luminosity Lx emitted 
from the surface of the neutron star is calculated with 
equation (5) by 

1.8 X 10^^- 



Lx = exmiAL 



Edd 



mi.4mi6 erg s 



(16) 

where ex = GMx/{Rxc^) ~ 0.2 represents the mass-to- 
energy conversion efficiency of the neutron star. By com- 
paring equation (15) with equation (16), we note that 
LxS^^ADAFj unless a is extremely small (a^ 10~^). This 
suggests that virtually all the observed X-ray luminosity 
arises from the polar cap regions of the neutron star. 

In this section, as a model for normal X-ray outbursts, 
we have numerically studied the mass-accretion rate via 
RIAFs of material captured from truncated Be disks, for 
which there is a gap between the disk and the neutron 
star orbit. No analytical approach was adopted here, be- 
cause mass comes from a tidally elongated part of the 
disk, where the density and velocity distribution can be 
obtained only numerically. In the next section, where we 
consider a mechanism for giant X-ray outbursts, we will 
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see that direct collision of the neutron star with a Be disk 
extending beyond the neutron star orbit can also cause 
smaller-scale outbursts with the luminosity at a level of 
normal outbursts. In studying such interaction, we will 
adopt an analytical approach as well as a numerical one 
for better understanding of physics at work there. 

4. Accretion Mechanism for Giant X-ray 
Outbursts 

In the previous section, we proposed that normal X- 
ray outbursts are caused by ADAFs of material trans- 
ferred from a truncated Be disk. The peak mass-transfer 
rate obtained by numerical simulations was in the range 
-|^gi6-i7 g g-i^ This rate, together with the accretion 
timescale of RIAFs, gives rise to the peak X-ray lumi- 
nosity of the order of lO'^^ergs"^. As summarized below, 
however, it is about one order of magnitude lower than 
that of giant X-ray outbursts. In other words, the mass- 
transfer rate comparable to 10^^~^*g s~^ is required to ex- 
plain the peak X-ray luminosity of giant X-ray outbursts. 

At the first glance, it seems very unlikely for a trun- 
cated Be disk to provide such a high mass-transfer rate, 
because there is a gap between the neutron star orbit and 
the Be disk. However, in misaligned systems, where the 
tidal torque is weaker than in coplanar systems, the trun- 
cation radius could be larger than the periastron separa- 
tion, depending on the tilt angle and the azimuth of tilt. 
If such a disk is warped towards the orbital plane and 
if its direction happens to be near the periastron direc- 
tion, the neutron star could capture gas at a high enough 
rate for a giant X-ray outburst, when it passages through 
the warped part of the Be disk. Indeed, disk warping 
episodes associated with giant X-ray outbursts have been 
observed, as mentioned in section 1 and also summarized 
in section 4.1. 

There are two possible mechanisms to make a Be 
disk warped; the radiation driven warping (Pringle 1996; 
Porter 1999) and the tidally driven warping (Martin et 
al. 2011). Unfortunately, our SPH code cannot currently 
handle either of these mechanisms, although there are 
SPH codes in which one of these mechanisms is imple- 
mented (Foulkes et al. 2010 for radiation-driven warping 
simulations; Lodato & Price 2010 for tidal warping sim- 
ulations). Hence, as a first step towards modeling giant 
X-ray outbursts, we study below the interaction between 
the neutron star and a rigidly tilted disk that extends 
beyond the neutron star orbit. After summarizing obser- 
vational characteristics of giant X-ray outbursts, we first 
consider the Bondi-Hoyle-Lyttleton accretion of Be-disk 
material as a simplest approach for modeling the passage 
of the neutron star through the Be disk, and then numer- 
ically study accretion fiows from differentially rotating Be 
disks for more detailed modeling. Whereas numerical sim- 
ulations enable us to quantitatively compare results, the 
analytical approach provides us qualitatively better un- 
derstanding of physics behind the numerical result. We 
also discuss the possibility that the accretion flows be- 
come supercritical if the Be disk density is very high. 



Observational Characteristics of Giant Outbursts 

The observed features of giant X-ray outbursts are 
summarized as follows: The X-ray luminosity increases 
by a factor ^10^ with respect to the quiescent state 
{Lx ^ 10'^^ crgs^^). It sometimes rises up by a factor 
of 10^~^ to a level comparable with the Eddington lumi- 
nosity of the neutron star. The duration time is several 
tens of days O.SPorb, sometimes over one orbital pe- 
riod). No orbital modulation has been detected. The 
fact that the large spin-up of the neutron star is seen 
in giant X-ray outbursts provides a strong evidence for 
the presence of a transient accretion disk during the giant 
outburst. Quasi-periodic oscillations (QPOs) detected in 
some systems (e.g., Takeshima et al. 1994; Finger et al. 
1996; Heindl et al. 1999) also support the transient accre- 
tion disk scenario. 

In some systems, the warping of Be disk associated with 
the giant outburst has been observed (Negueruela et al. 
2001; Reig et al. 2007; Moritani et al. 2011). Negueruela 
et al. (2001) and Reig et al. (2007) has made long-term 
monitoring observations of the Ha emission line profile 
from the optical counterpart of 4U 0115-1-634 in 1995- 
2005, and found that before and during giant X-ray out- 
bursts, it changed from usual double-peaked profile to a 
single-peaked or shell-line profile on a timescale of a year 
or so. As mentioned in section 1, this profile variability is 
most likely to originate from the precession of a warped Be 
disk. Moreover, high-resolution Ha profiles of the optical 
counterpart of A 0535+262 observed during the giant X- 
ray outburst in 2009 had an enhanced component, which 
together with the double-peaked component, appeared as 
a triple peak. Further analysis of these profiles indicates 
that the Be disk became tidally warped (Martin et al. 
2011) and intersected with the binary orbital plane near 
the periastron during the giant X-ray outburst (Moritani 
ct al. 2012, in prep). 

4.2. Bondi-Hoyle-Lyttleton Accretion 

We consider the Bondi-Hoyle-Lyttleton (BHL) accre- 
tion process where the neutron star captures gas from an 
initially uniform flow of density p and relative speed Urei- 
The BHL accretion was first analyzed by Bondi & Hoyle 
(1944) and have been extensively studied since then (see 
Edgar 2004 for references and a modern review). 

In the classical BHL accretion, the material with initial 
impact parameter smaller than a critical value 

, 2GMx 1, / v,^i 

b= ^ «3.7x 10i2mi4 — r cm 17 

■ V 100 km s"V 

has negative total energy, and is eventually accreted on to 
the central object. This leads to a simple BHL formula 
for the expected mass accretion rate, 

= 4.4 X 10^ V? 4 I -r^ ) 
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Fig. 6. Schematic diagram of a misaligned Bc/X-ray binary. 

In order to emulate the accretion process from an ex- 
tended tilted disk, we adopt the disk gas density and veloc- 
ity evaluated at the position of the neutron star, assuming 
the Be disk to cross the orbit at this position (sec figure 6 
for the model geometry). In the calculation, the Be disk 
is assumed to rotate at the Kcplcrian speed and have the 
density distribution in the form, 



P'- 



X 1 — 



(19) 



^^g cm ^ and 



with typical Be disk parameters, po = 10 
n = 7/2. Here, r is the distance from the Be star. For 
simplicity, we neglect the velocity shear and density gra- 
dient around the neutron star. 

As an example, in the top panels of figure 7, we show 
the relative speed Wroi between the neutron star and a Be 
disk with tilt angle /3 = 45° as a function of azimuth of the 
ascending node 7-|-7r/2. The left panel is for A 0535-1-262, 
while the right panel is for 4U 0115-1-634. In both calcu- 
lations, we assumed that the disk rotates in the prograde 
direction. Note that the plot shows the dependence of 
the relative velocity on the azimuth of tilt, 7, of the Be 
disk, not on the orbital phase of the neutron star. For 
7-|-7r/2 = 0, the neutron star passes through the Be disk 
at the apastron, while for 7 + 7r/2 = tt the periastron is 
the location of the disk transit. 

We mentioned above that the material with initial im- 
pact parameter smaller than b accretes to the neutron 
star. In binaries, however, only material inside the Roche 
radius of the neutron star, tl, given by equation (4), can 
be accreted. Thus, in the BHL accretion in binaries, the 
accretion rate is given by 



BHL ■ 



Mo X min 



(20) 



where Mq is the classical BHL accretion rate given by 
equation (18). The middle panels of figure 7 compare 
these two radii, b (blue line) and tl (red line). The solid 
part of these lines indicates the effective accretion radius, 
i.e., the smaller radius of h and tl, inside which the gas 
accretes on to the neutron star. We note that for /3 = 45°, 
the accretion radius in A 0535-1-262 is controlled by b when 
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Fig. 8. Mass-accretion rate as a function of tilt angle f). 
The solid and the dashed lines denote analytical Mbhl 
for A 0535-1-262 and 4U 0115-1-634, respectively, while the 
red and the blue circles denote the peak mass-accretion 
rate obtained by numerical simulations of A 0535-1-262 and 
4U 0115-1-634, respectively. In these calculations, the az- 
imuth of tilt, 7, is fixed at 90°, for which the neu- 
tron star passes through the Be disk at the periastron. 

the neutron star passes through the disk around the pe- 
riastron, while it is equal to tl otherwise. In our model, 
this transition occurs only for moderate values of tilt an- 
gle. This is also the case in 4U 0115-1-634, although in this 
smaller system with higher relative velocity, the accretion 
radius in the particular case of /? = 45° is controlled by b 
for any value of azimuth of tilt 7. For /3 -C 45°, because 
of small relative velocity, the critical radius b for the BHL 
accretion exceeds the Roche radius tl at any azimuth of 
the disk plane, so that the accretion radius is always given 
by tl. On the other hand, for [5 » 45° (including all cases 
of retrograde disk rotation), the relative velocity is so high 
that the BHL accretion radius b is always smaller than the 
Roche radius tl, as in the case of stellar wind accretion. 

The bottom panels of figure 7 present the accretion rate 
.^^BHL calculated by using equation (20). Note that in 
A 0535-1-262, the mass-accretion rate of ^ 10^ g s^^, and 
hence the accretion luminosity of ^ 10'^^ ergs~^, is ob- 
tained for a wide range of azimuth of the misaligned disk. 
In contrast, in 4U 0115-1-634, such a high accretion rate is 
obtained only if the Be disk density is significantly higher 
than the typical density. 

Figure 8 shows the mass accretion rate as a func- 
tion of tilt angle (3. In these calculations, the azimuth 
of tilt 7 was fixed at 90°, for which the neutron star 
passes through the Be disk at the periastron. The solid 
and the dashed lines denote Mrhl for A 0535-1-262 and 
4U 0115-1-634, respectively. From the figure, we note that 
in A 0535-1-262 Mbhl ^ 10'^^ erg s~^ is obtained for a very 
wide range of tilt angle for typical Be disk parameters 
(po 10-"gcm-3 and n = 7/2), while in 4U 0115-^634 
it is reached only for a small range of /3 as long as the 
typical Be disk parameters are adopted. This difference 
arises from the fact that larger systems (A 0535-f262 in 
our sample systems) have larger Roche radii, which make 
the accretion rate higher. From figure 8, we also note that 





I ' I I ' I 

-0.5 0.5 1 1.5 -0.5 0.5 1 1.5 

(y + nl2)l 2% (j + nl2)l 2n 

Fig. 7. Characteristic quantities of the BHL accretion as a function of ascending-node azimuth 7 + 77/2, at which the neutron star 
passes through the Be disk: (a) A 0535+262 and (b) 4U 0115+634. Top panels: Relative velocity between the neutron star and 
the Be disk. Middle panels: Comparison between BHL accretion radius b (dotted line) and Roche radius rL (dashed line) of the 
neutron star. Bottom panels: Mass-accretion rate Mbhl given by equation (20). In these calculations, the tilt angle /3 is fixed 
at 45°. The base density and the density gradient index of the Be disk are assumed to be 10~^^gcm~'' and 7/2, respectively. 



the highest accretion rate occurs at /3 '~ 30 — 40°. This an- 
gle is shghtly smaller in 4U 0115-^634 than in A 0535-^262. 
This is because small systems have relatively larger rela- 
tive velocity. 

In this subsection, we have considered the BHL accre- 
tion as a model for the direct accretion from the mis- 
aligned Be disk. However, the Be disk is differentially ro- 
tating, so that accretion occurs in a sheared flow. The 
finite thickness of the disk and the eccentric orbit of 
the neutron star should also affect the accretion rate. 
Therefore, the BHL accretion assuming a uniform back- 
ground flow should be taken as a first order model. In 
the next section, we numerically examine accretion flows 
in differentially rotating Be disks, using SPH simulations. 



4-.3. Accretion from Differentially Rotating Be Disks 

Although very simplified, the above analytical model 
provides a useful information about the mass captured 
by the neutron star when it passes through a tilted Be 
disk. The resultant mass-capture rate is of the order of 
lO^^ergs"^ for a wide range of tilt angles. This rate is 
comparable to the mass accretion rate required for the 
observed luminosity of giant X-ray outbursts. 

As seen in section 2.2.2, however, high mass-capture 
rates do not necessarily mean high mass-accretion rates. 
If the accretion time is comparable with or longer than 
the orbital period, the accretion luminosity will stay more 
or less constant at a level much lower than that calculated 
from the peak mass-capture rate. Hence, it is important to 
examine the accretion timescale of mass captured by the 
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a) A 0535+262: /3 = 20° (b) A 0535+262: /3 = 40° 



neutron star. It is, however, not a trivial task to estimate 
the accretion rate and timescale, given the presence of 
differential rotation in the Be disk and the eccentricity of 
the neutron star orbit. For quantitative study of these 
quantities, we carried out SPH simulations of accretion 
processes from misaligned Be disks for a limited range of 
parameter space for A 0535+262 and 4U 0115+634. 
4-3.1. Numerical Model 

The initial configuration of simulations are as follows. 
For simplicity, we consider a rigidly-tilted, isothermal Be 
disk to model the interaction between the neutron star 
and a misaligned disk. In order to simulate direct colli- 
sion of the neutron star with the Be disk with good spa- 
tial resolution, we concentrated SPH particles in a ring- 
like region across the periastron separation r = (1 — e)a, 
from r ^ 0.65(1 — e)a to r ^ 1.2(1 — e)a. We distributed 
SPH particles such that the Be disk has the equatorial 
density distribution in the form of equation (19) with 
Po = 10~^^g cm~'^ and n = 7/2 and is in hydrostatic equi- 
librium in the direction normal to the disk mid-plane. The 
disk temperature is set to 0.6reff as in section 2.2.1. We 
tilted the Be disk about the semi-major axis (i.e., 7 = 90°) 
to study the strongest interaction case for each tilt an- 
gle. The accretion radius of the neutron star is set to 
O.Ola « 2 X 10^i?x- The initial number of SPH particles is 
7 X 10^. We have run each simulation from O.lPorb prior 
to the periastron through O.lPoib after it. 
4.3.2. Numerical Results 

We have run simulations for five different tilt angles 
/? = 0°, 20°, 40°, 60°, and 90° for each system. The peak 
accretion rates obtained in these simulations are plotted 
by filled circles in Figure 8. In contrast to the analytical 
result, the simulated peak accretion rate monotonically 
decreases with increasing misaligned angle. Compared 
with the analytical rate, the simulated rate is significantly 
high in the coplanar (/3 = 0) case, comparable for /3 = 20° 
and 90°, and significantly low for /3 = 40° and 60°. 

In order to understand the cause(s) of these differences, 
we compare simulation results for tilt angles (3 = 20° and 
P = 40° in more detail in figures 9 (A 0535+262) and 
10 (4U 0115+634). In each figure, upper panels show 
the orbital-phase dependence of the mass accretion rate 
on to the neutron star. In each panel, the thick (blue) 
line denotes the rate of mass inflow through the inner 
simulation boundary (hereafter, the simulated accretion 
rate). 

For /3 = 20°, the simulated accretion rate rises rapidly, 
has a pre-periastron peak, and then decays slowly. On 
the other hand, for /3 = 40°, there are two peaks: The 
lower pre-periastron peak is followed by the higher post- 
periastron peak. After the post-periastron peak, the rate 
decreases gradually. In both cases, the initial accretion is 
of the BHL type, but later (after the neutron star passes 
through the disk) the accretion is via an accretion disk 
formed by material that catches up with the neutron star 
from behind, which has significantly larger specific angular 
momentum than the martial accreted before and during 
the disk transit. The initial BHL accretion and the sub- 
sequent accretion via the accretion disk form the rapidly 
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Fig. 9. Simulation result for tiie BHL accretion from a tilted 
Be disk in A 0535+262: (a) /3 = 20° and (b) ^ = 40°, where 
/3 is the tilt angle. The azimuth of tilt 7 is fixed at 90°, 
for which the neutron star passes through the Be disk at 
the periastron. Each simulation has run from orbital phase 
—0.1 to -1-0. 1. Upper panels: The orbital-phase dependence of 
the mass-inflow rate through the inner simulation boundary 
[thick (blue) lines] and the estimated accretion rate at the neu- 
tron star radius [thin (red) lines)]. The right axis shows the 
X-ray luminosity corresponding to these rates. Lower pan- 
els: The orbital-phase dependence of the accretion timescale 
[thick (red) lines]. The standard and the ADAF accretion 
timescales are also shown by the upper and the lower thin 
(black) lines, respectively. The horizontal (red) dashed line 
in each lower panel denotes the mass- weighted average of the 
accretion timescale. 

(a) 4U 0115+634: /3 = 20° (b) 4U 0115+634: /3 = 40° 
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Fig. 10. Simulation result for the BHL accretion from a 
tilted Be disk in 4U 0115+634: (a) ^ = 20° and (b) ^ = 40°. 
The azimuth of tilt 7 is fixed at 90°, for which the neutron 
star passes through the Be disk at the periastron. Each sim- 
ulation has run from orbital phase —0.1 to +0.1. The format 
of the figure is the same as that of figure 9. 

rising part and the slowly decline part of the simulated 
accretion rate profile, respectively. 

As seen in figure 8, the analytical model predicts a 
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higher accretion rate for /3 ^ 30° - 40° than for /3 = 20°. 
The simulations, however, have revealed the opposite 
trend. The simulated accretion rate for (3 = 20° is 2- 
3 times as high as that for /? = 40° in both systems. 
Examining the interaction more closely, we have found 
that in the case of small tilt angles, the tidal torque by 
the neutron star starts exciting a tidal stream in the Be 
disk significantly before the disk crossing event, which col- 
lides with the neutron star at a very small impact pa- 
rameter. Since the density in the tidal stream is much 
higher than in the background disk gas, this collision en- 
hances the mass-accretion rate significantly. In contrast, 
for (3 ^ 40°, where the tidal torque is weaker, the tidal 
stream grows more slowly. By the time it becomes strong, 
the neutron star is already passing the periastron and the 
streaming gas, which has a relatively large impact param- 
eter, forms an accretion disk, not directly running on to 
the neutron star, thus unable to significantly enhance the 
mass-accretion rate. 

Since the accretion radius set in these simulations is 
O.Ola ^ 2 X 10^i?x and still way far from the neutron star, 
the accretion time the flow takes from this radius to the 
neutron star surface could significantly modify the ampli- 
tude and shape of the accretion rate profile. Thus, we 
estimated the accretion rate at the neutron star radius 
(hereafter, the final accretion rate) using the same pro- 
cedure described in section 2.2.2. In deriving the final 
accretion rate, we assumed that the accretion flow is of 
standard-disk type in the period of m > 1 (or mig > 20) 
and afterwards, while it is of ADAF type otherwise. In 
the former (latter) case, the accretion takes place in the 
viscous (ADAF) timescale given by equation (3) [equa- 
tion (13)], with the circularization radius derived from 
the specific angular momentum of accreted SPH parti- 
cles. The phase dependence of the accretion timescale, 
iacc, is shown by the thick (red) lines in the lower pan- 
els of figures 9 and 10, where for comparison purpose, 
the viscous and ADAF timescales are also denoted by the 
upper and lower thin lines, respectively. The resultant, 
final accretion rate is shown by the thin (red) lines in 
the upper panels. From the figures we note that the fi- 
nal accretion rate profile changes little from the simulated 
profile at r = O.Ola, except for cases where the accretion 
rate exceeds the upper limit for ADAF solutions, so that 
the accretion fiow is of standard-disk type. For typical 
Be disk parameters, the latter case occurs only in (3 < 20° 
simulations for A 0535-1-262. 

The above simulations of accretion from differentially 
rotating Be disks have important observational implica- 
tions. The peak accretion rates obtained for /? ^ 20° and 
typical Be disk density {p^n ^ 1) agree with those re- 
quired for the observed X-ray luminosity of giant X-ray 
outbursts. In other words. Be disks of typical density have 
potential to cause giant outbursts if the tilt angle is small. 
On the other hand, if the tilt angle is large, say (3 ^ 40°, 
the Be disk has to be several times denser than typical 
in order to supply large enough amount of gas for giant 
X-ray outbursts. Otherwise, the outbursts will appear as 
normal X-ray outbursts. 



4-4- Supercritical Accretion 

In the previous section, we assumed that the accretion 
of mass captured by the BHL-type accretion occurs via the 
standard disk or optically thin ADAF. This assumption is 
reasonable as long as the mass-capture rate is of the order 
of 10^''g s~^ or less. If the mass-capture rate exceeds the 
Eddington rate (~ lO^^g s~^), however, the accretion flow 
drastically changes to being supercritical. In this section, 
we investigate the possibility that Bc/X-ray binaries have 
supercritical accretion flows. 

Recent X-ray observations have revealed the existence 
of bright X-ray sources such as GRS 1915-1-105 in our 
galaxy with luminosities over the Eddington luminosity 
(c.f.. Done et al. 2007). Such large luminosities can be 
explained by a supercritical accretion flow on to a black 
hole or neutron star (c.f., Watarai et al. 2000). 

In an optically thick accretion flow with an accretion 
rate much higher than the Eddington accretion rate, the 
photons are trapped and restored as the entropy in the 
accreting gas without being radiated away. Such a flow is 
called a slim disk, where the advective cooling dominates 
the radiative cooling. The slim disk model was first pro- 
posed by Abramowicz et al. (1988) as a thermally stable, 
advection-dominated, supercritical accretion flow. Then, 
its self-similar solution has been derived by Wang & Zhou 
(1999) and Watarai (2006). 

In what follows, we flrst describe the condition that 
the supercritical accretion occurs in the accretion flow 
around the neutron star in the context of Be/X-ray bina- 
ries. We then derive the accretion timescale based on the 
self-similar solutions of slim disk model (Watarai 2006). 
We flnally evaluate the Bolometric X-ray luminosity emit- 
ted from the slim disk. 

4.4-1- Condition for Supercritical Accretion 

At luminosities close to or higher than the Eddington 
luminosity, the scale-height of the accretion disk is compa- 
rable to the radius, i.e., H/r ^ 1. Hence, in the following 
analysis, we do not distinguish between spherically sym- 
metric accretion and disk accretion. 

In a very optically thick medium, photons are trapped 
if their diffusion timescale ^diff is longer than the accretion 
timescale iacc- If the medium is disk-like, tdiff of photons 
near the equatorial plane is given by 



idiff 



MP 



3Ht 3HkY. 



2c 



(21) 



(e.g., Kato et al. 2008), where A is the mean free path of 
photons, T is the vertical optical depth, k is the opacity, 
and E is the surface density. In deriving the last equation, 
we used r ~ kE/2. As the opacity k, we use that of the 
electron scattering, Kos (^0.4). 

On the other hand, the accretion timescale is written as 



t. 



Vr 



27rr2S 
M 



(22) 



Equating equation (21) with equation (22) and using 
equation (5), we have the photon-trapping radius rtrap 
as 
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3 H Mc^ 

nrap = 7; — 7 '^s « 3.2 X 10 cm, (23) 

where H /r ~ 1 at high luminosities. In order for the su- 
percritical accretion to occur, rtrap should be larger than 
the magnetospheric radius tm, from which we have the 
following condition: 

mi6 > mi6,sc, (24) 

where we define 

mi6,,e = 1.7 X l{y'k'/''m-'-J''^l%\ (25) 

From equation (25), the mass-capture rate must be higher 
than ~ 2 X lO-'^^g s~^ in order that the accretion flow be- 
comes supercritical before it reaches the magnetosphere 
of the neutron star. To realize such a huge mass- 
capture rate, the base density of the Be disk should 
be by two orders of magnitude higher than the typical 
Po = 10"^^gcm~'^. Note that this is not an unrealistic 
condition. For instance, figure 9(a) shows that in the case 
of A 0535-1-262 with 13 = 20° and 7 = 90°, the accretion 
flow can be supercritical at orbital phases from —0.04 to 
0.02 with respect to the periastron, if the base density is 
10~^gcm~^ or p_ii = 100. 
4-. 4-2. Observational Implications 

In this subsection, we evaluate the Bolometric X-ray 
luminosity emitted from the slim disk. We assume that 
the condition to > mgc is satisfied so that the slim disk 
state holds from the photon-trapping radius down to the 
magneto-spheric radius. 

According to a self-similar solution of the slim disk 
(Watarai 2006), the radial velocity is given by 

«3.2 X 10Vao.i^'"^/^cms"\ (26) 

where we recall that / = Qadv/Q^- Then, the accretion 
timescale is calculated as 

tscAF = ^ « 2.7 X 10- V«o:i™^y'^?6' s- (27) 

The distribution of effective temperature of the slim disk 
is given by the self-similar solution as 

reff«4.9xl0V'/V-yVi/2K (28) 

(Watarai 2006), so that the cooling rate of the slim disk 
Qrad is given by 

Qiad = 2(7 T^^ 

w 1.1 x IQ^^f^/'^m^lf-'^ ergs-i cm'^. (29) 

The Bolometric X-ray luminosity emitted from the slim 
disk is then estimated as 

isliin = 27rrg / rQraddf 

J f-M 

« 1.2 X lO^'f/^mi 4ln ( ^) eTgs~\ (30) 

From equation (27), the timescale of supercritical accre- 
tion is much shorter than the orbital period of any Be/X- 
ray binaries. Hence, if the mass-capture rate satisfies the 
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condition given by equation (24), the resultant supercrit- 
ical accretion fiow will have the same final accretion rate 
profile as those shown by blue lines in the upper panels of 
figures 9 and 10 with the amplitude proportional to the 
base density of the Be disk. From equation (30), the lumi- 
nosity is then estimated to be of the order of lO^^erg 
even if rtrap/ru ~ 10. Note that the light curves of super- 
critical accretion flows in Be/X-ray binaries are thus quite 
different from those expected for standard disk accretion. 

5. Discussion 

In Be/X-ray binaries, the Be disk is tidally /resonantly 
truncated by the gravitational force of the neutron star, 
except in highly eccentric systems and those with Be 
disks highly inclined from the orbital plane (Okazaki & 
Negueruela 2001). Such a truncation produces a gap be- 
tween the Be disk and the neutron star orbit. In a system 
with a gap, the neutron star captures gas from a tidally 
elongated part of the disk. Since the tidal deformation 
is a dynamical process and there occurs no redistribution 
of the angular momentum by viscosity, the rotation ve- 
locity in the elongated part decreases more rapidly with 
radius than in a Keplerian disk. This makes the rela- 
tive velocity between the neutron star and the captured 
material larger than in the case where the neutron star 
captures mass from an extended Keplerian disk without a 
gap. As a result, the material captured from a truncated 
Be disk forms a large accretion disk around the neutron 
star. If the accretion disk is radiative-cooling dominated, 
its accretion timescale easily exceeds the orbital period. 
In such a situation the system will exhibit no rapid nor 
large X-ray flux changes seen in normal X-ray outbursts. 
In fact, in section 2, we have shown that this is the case 
for the Be/X-ray binaries A 0535-^262 and 4U 01154-634, 
by using results from numerical simulations. 

One might argue that the magnetosphere of the neu- 
tron star controls accretion flows, causing rapid and large 
variation of X-ray luminosity. When the magnetospheric 
radius is larger than the coronation radius of the accreting 
matter, the matter is expelled by the rotating magneto- 
sphere (the propeller effect) and the systems is in the qui- 
escent state. The gate opens only if the Be disk supplies 
large enough mass so as to the resultant accretion flow 
can push back the magnetosphere inside the corotation 
radius. In fact, in 4U 0115+63, a sudden transition from 
the quiescence to a normal outburst has been witnessed, 
which is consistent with this idea (Campana et al. 2001). 

However, as shown in section 2, A 0535-1-262 is almost 
always in the direct accretion regime (flgure 3), because 
the propeller effect does not work for this slowly rotat- 
ing X-ray pulsar. The system is therefore one of best 
candidates to study more direct relationship between the 
accretion flows on to the neutron star and the resultant 
X-ray variability. Given that standard disks have accre- 
tion timescales much longer than the orbital period, we 
can safely rule out the possibility that the material from 
the truncated Be disk forms the standard disk around the 
neutron star. 
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The problem with rapid and large X-ray variability 
is resolved if the accretion flow is radiatively inefficient 
(RIAF). In section 3, using results from numerical simu- 
lations, we have shown that mass transfer from truncated 
Be disks with typical density parameters gives rise to ac- 
cretion rates in agreement with those expected from X-ray 
luminosity of normal outbursts, and that these rates are 
significantly lower than the critical rate above which no 
RIAF solution exists. The accretion timescale is also con- 
sistent with the duration of normal outbursts. 

The RIAF model for normal X-ray outbursts opens a 
new research field on Be/X-ray binaries. In this paper, we 
assumed that the magnetospheric radius for self-similar 
ADAF solutions is about the same as that for disk ac- 
cretion. Given the magnetospheric interaction affects the 
accretion on to the neutron star in many ways, including 
even the possibility of launching jets, it is very important 
to investigate it in detail. There is also a possibility to de- 
tect high energy emission from RIAFs, e.g., synchrotron 
emission in X-rays from accelerated particles in the layer 
above the magnetosphere and the GeV/TeV emission via 
the inverse Compton scattering of soft photons from the 
Be star and/or the accretion disk, although Fermi and 
VERITAS have detected no emission above 0.1 GeV from 
A 0535-f 262 during the 2009 giant outburst (Acciari et al. 
2011). We will study it in a forthcoming paper. 

Another interesting possibility also emerges from the 
RIAF model for normal X-ray outbursts. In section 3, we 
have used the numerical results for a typical Be disk with 
base density of 10~^^gcm~'^ and density gradient index 
n = 7/2. There is, however, a large range of observed 
base density (10~^^g cm^'^ — 10^^°g cm^'^) and density 
gradient index (2.5-4.5) (e.g., Jones et al. 2008). If the 
mass transfer rate from the Be disk exceeds the critical 
rate above which no RIAF solution exists, a standard ac- 
cretion disk would form. Thus, the transition between 
standard disk state and RIAF state can occur, depending 
on the base density and density gradient of the Be disk. 
Obscrvationally, this transition will be seen as a sudden 
drop/rise of X-ray luminosity by about one order of mag- 
nitude. Since the mass-transfer rate decreases with in- 
creasing angle between the Be disk and the binary orbital 
plane, the coplanar Be-disk systems are most favorable to 
see this transition. 

In section 4, we have proposed that giant X-ray out- 
bursts can be triggered by the BHL accretion of material 
from a tilted Be disk. Although the RIAF model with 
truncated Be disks works for normal X-ray outbursts, it 
cannot be applied to giant X-ray outbursts by two rea- 
sons. One reason is that because of the disk truncation, 
the mass-transfer rate is lower than that required for gi- 
ant X-ray outbursts, unless the Be disk density is much 
higher or its distribution is much flatter than typical val- 
ues. The other, more crucial reason is that the mass- 
accretion rate estimated from the luminosity of bright gi- 
ant X-ray outbursts exceeds by a factor of few the highest 
mass-accretion rate for RIAFs. Moreover, there is ob- 
servational evidence that the Be disk is warped when gi- 
ant outbursts occur (Negueruela et al. 2001; Reig et al. 



2007; Moritani et al. 2011). 

Our model for giant X-ray outbursts is based on the idea 
that the warped shape of the Be disk naturally results in 
the enhanced mass accretion, if the disk is inclined from 
the orbital plane and its warped part gets across the orbit 
of the neutron star from above or below the orbital plane. 
In real systems, the accretion rate should depend on the 
local density and velocity distribution of gas around the 
neutron star. For a warped disk picture, however, no such 
information is easily available. Therefore, as a first order 
of approximation, we have modeled the warped Be disk 
by a rigidly-tilted, Keplerian disk with a power-law den- 
sity distribution in the radial direction. The Be disk was 
assumed to have an outer radius larger than the perias- 
tron separation. With a simplified analytical treatment, 
we found that the accretion radius of the neutron star 
is limited by the Roche lobe radius for small tilt angles, 
while for large tilt angles the neutron star accretes mate- 
rial inside the BHL accretion radius. 

Numerical simulations also revealed a new feature of 
the BHL accretion in Be/X-ray binaries. For Be disks 
with small tilt angles {(3 < 20°), the tidal torque of the 
neutron star rapidly excites a coUimated gas flow in the Be 
disk, and the neutron star collides with this tidal stream 
at a very small impact parameter. Since the density in 
the tidal stream is much higher than in the background 
disk gas, this collision greatly enhances the mass-accretion 
rate. This does not happen for larger tilt angles (/3 ^ 40°). 
Consequently, the accretion rate is much higher for small 
tilt angles than for large tilt angles. 

Using simulation data, we have estimated the time- 
dependent accretion rate at the neutron star radius. It 
turns out that for typical Be disk parameters, only long 
period systems, such as A 0535-1-262, with small tilt angles 
{f3 ^ 20°) can have mass-accretion rates as high as those 
estimated for giant X-ray outbursts. In order for systems 
with large tilt angles {(3 ^40°) and/or short periods such 
as 4U 0115-1-634 to exhibit giant X-ray outbursts, the Be 
disk density has to be higher than the typical one by a 
factor of several or more around the neutron star orbit. 
Otherwise, the outburst luminosity of these systems falls 
into the range of normal X-ray outbursts. 

In this paper, we have assumed the upper limit for 
the ADAF solutions to be to = M/(LEdd/c^) = 1 (or 
A/ 2 • 10^^ g s~^). With this criterion, the accretion 
flow is radiative-cooling dominated (i.e., a standard disk) 
only for relatively bright (^ 4- 10^^ erg s""'^) giant out- 
bursts. For less bright (^ 4- lO^^ergs"^) giant outbursts 
and all normal outbursts, the accretion flow is advection- 
dominated. The type of the accretion flow, standard-disk 
type or ADAF type, thus depends only on whether the 
mass-accretion rate is higher than the limiting rate for 
ADAF solutions or not. 

From the sensitiveness of the accretion rate on the lo- 
cal tilt angle, a working hypothesis comes up for a series 
of X-ray activities consisting of a few normal X-ray out- 
bursts and a giant X-ray outburst. This series of X-ray 
activity has sometimes been observed in Be/X-ray bina- 
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ries. As mentioned in section 4.1, 4U 0115+634 has shown 
observational evidence of precession of a warped Be disk 
before and during giant X-ray outburst (Ncgueruela et al. 
2001; Reig et al. 2007). When a warped disk precesses, it 
is likely that the local tilt angle, i.e., the angle between 
the neutron star's motion and the motion of material in 
the Be disk, varies from cycle to cycle. If the tilt angle 
is large at the beginning, the system exhibits an X-ray 
outburst at a level considered to be a normal outburst. 
If the warped structure survives over several orbital cy- 
cles, there will be a good chance for the neutron star to 
pass through the disk at a small tilt angle in due course. 
This encounter will give rise to a giant X-ray outburst. 
Later on, the system will show a few outbursts at the level 
of normal X-ray outbursts until the azimuth of the disk 
plane (7-|-7r/2) moves too far away from the periastron. 
In this way, our model for giant X-ray outbursts natu- 
rally, though qualitatively, explains a sequence of X-ray 
outbursts, i.e., the occurrence of a few normal X-ray out- 
bursts before and after a giant X-ray outburst, which has 
been observed in 4U 0115-^634, A 0535-^262, and some 
other Be/X-ray binaries. 

Finally, there is possibility of supercritical accre- 
tion in Bc/X-ray binaries, as discussed in section 4.4. 
Throughout the paper, we fixed the base density and 
the density gradient index of the Be disk at 10""'^^ g cm~'^ 
and 7/2, respectively. As mentioned previously, however, 
there is a large scatter in the distribution of the Be disk 
base density and its density gradient. For Be disks with 
very high base density and/or significantly fiatter den- 
sity distribution, the model predicts that the accretion 
becomes supercritical, triggering super-Eddington lumi- 
nosity. As long as we know, no Be/X-ray binaries have 
shown super-Eddington luminosity. Nevertheless, we see 
no mechanism that can prevent Be disks from supplying 
mass at supercritical rates. If the supercritical accretion 
takes place, the luminosity will exceed or be comparable to 
the Eddington luminosity and photon trapping will veil all 
features related to magnetic fields of the neutron star such 
as X-ray pulsation and cyclotron resonance scattering fea- 
tures. Once such characteristic changes arc observed, a 
fertile research field will arise from it. It is therefore im- 
portant to carry out more detailed modeling of supercrit- 
ical accretion in Be/X-ray binaries. Be/X-ray binaries 
thus provide a unique opportunity where in one system 
one can study physics of various types of accretion, in- 
cluding standard disks, Bondi-Hoylc-Lyttlcton accretion, 
RIAFs, and even supercritical accretion fiows, in relation 
to the geometry and dynamics of the mass donor. 

6. Conclusions 

We have studied the origins of normal and giant X-ray 
outbursts in Be/X-ray binaries. Based on the results from 
3D SPH simulations, we first showed that the accretion via 
the standard disk cannot produce the observed rapid and 
large X-ray fiux variation. Then, by using analytical stud- 
ies as well as numerical ones, we proposed the following 
scenarios for these two types of X-ray outbursts: 



1. The normal X-ray outbursts are caused by RIAFs 
from tidally truncated Be disks. The resultant lumi- 
nosity and time variability are consistent with those 
observed, because RIAFs have accretion timcscale 
much shorter than the orbital period, so that the 
rapid and large variation in the mass-transfer rate 
from the truncated Be disk is conserved during the 
accretion process. 

2. The giant X-ray outbursts occur in systems where 
the Be disk is misaligned with the binary orbital 
plane and is sufficiently developed. When the out- 
ermost part of such a disk is warped and crosses 
the orbit of the neutron star near the periastron, 
the neutron star can capture a large amount of gas 
via the BHL accretion. This process results in the 
formation of a very small standard disk or ADAF, 
depending on the mass-capture rate, in which the 
accretion timescale is much shorter than the orbital 
period. 

In this paper, we focused on the mass-transfer process 
from the Be disk to the neutron star and the subsequent 
accretion processes. There are, however, observational 
hints that the giant X-ray outburst is not a single inde- 
pendent event, but is a highlight of a much longer cycle of 
events. For instance, a giant X-ray outburst is often ac- 
companied by smaller scale X-ray outbursts, and the Be 
disk emission starts declining sometime before the giant 
outburst. In the next paper, we will study how the giant 
X-ray outburst and these associated events are related to 
the evolutionary cycle of the Be disk. 
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